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Effects  of  doping  level  and  compensation 
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ABSTRACT 

It  is  shown  that  the  thermal  ionization  energy  of  Mg  acceptors  in  GaN,  as  determined  by  temperature  dependent  Hall  effect 
measurements,  exhibits  the  usual  dependence  on  the  concentration  of  ionized  impurities,  as  seen  in  many  other 
semiconductors.  The  observed  difference  in  the  thermal  and  optical  ionization  energies  of  Mg  acceptors  can  be 
quantitatively  understood  based  on  a simple  electrostatic  interaction  model. 
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1.  INTRODUCTION 

Gallium  nitride  (GaN)  is  a promising  direct  wide  band  gap  semiconductor  for  potential  applications  in  visible  and  near- 
UV  optoelectronics  and  in  high  temperature  electronic  devices  (for  recent  reviews  see  e.g.  references1^).  In  fact  the 
invention  of  blue  semiconductor  lasers  made  from  gallium  nitride  has  been  one  of  the  key  technological  breakthroughs  of 
the  decade.  One  of  the  key  steps  forward  in  the  long  way  to  the  realization  of  these  promises  was  the  achievement  of 
controlled  p-type  doping  using  Mg  as  an  acceptor  by  Amano  et  al.5  and  later  by  Nakamura  et  al.6  One  of  the  basic 
characteristics  of  an  acceptor  dopant  is  its  ionization  energy.  The  ionization  energy  of  the  Mg  acceptor  as  determined  by 
optical  (low  temperature  photoluminescence)  methods  is  224±4  meV7  However  reported  values  of  thermal  ionization 
energies  for  the  Mg  acceptor  are  considerably  lower,  with  values  scattered  in  the  range  from  125  meV  to  170  meV  based 
on  temperature  dependent  Hall  effect  measurements  (see  e.g.  reference^),  and  about  155  to  165  meV  based  on  the 
temperature  dependence  of  the  intensity  of  donor-acceptor  pair  emission.^  These  large  differences  are  usually  ascribed  to 
the  ubiquitous  concentration  dependence  of  the  thermal  ionization  energy  of  impurity  centres  in  semiconductors,  however, 
up  to  now  a systematic  quantitative  analysis  has  not  yet  been  attempted.  For  device  applications,  however,  the  thermal 
ionization  energy  of  the  dopant  is  the  relevant  parameter,  because  it  controls  the  degree  of  ionization  of  the  dopant  centres, 
and  in  this  way  it  determines  the  available  free  charge  carrier  density,  dopant  efficiency,  etc.  Because  of  this  it  is  important 
to  understand  these  significant  differences. 

Here  I present  a quantitative  analysis  and  interpretation  of  the  observed  differences  in  the  thermal  and  optical 
ionization  energies  of  Mg  acceptors  in  GaN  on  the  basis  of  the  available  models  for  the  concentration  and  compensation 
degree  dependence  of  thermal  ionization  energies  of  impurity  centres.  It  is  proposed  that  these  differences  can  be 
satisfactorily  described  using  a simple  electrostatic  interaction  model  originally  due  to  Pearson  and  Bardeen  ^ and  Debye 
and  Conwell, 1 * and  later  discussed  and  generalized  by  Monecke  et  al.,1^  and  by  the  present  author.  13-15 


2.  MODELS  FOR  THE  CONCENTRATION  DEPENDENCE  OF  THE  THERMAL 

IONIZATION  ENERGY 

Since  the  pioneering  work  of  Pearson  and  Bardeen111  it  has  been  well  known  that  the  thermal  activation  energy  of  the 
impurities  decreases  with  increasing  concentration  of  the  impurity  centres.  Various  models1 0-12, 15  jea(j  t0  following 
generic  expression  for  the  dependence  of  the  thermal  ionization  energy  on  the  majority  impurity  concentration  (acceptors 
in  p-type  materials,  to  be  definitive)  and  on  the  compensation  degree: 
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EA  = EA(NA,K)  = EAo-f(K) 


(1) 


e2 

NA1/3 

47cej.£0 

here  EAo  is  the  ionization  energy  of  acceptors  at  infinite  dilution,  NA  is  the  acceptor  concentration,  K is  the  compensation 
degree  (K  = NjyNA),  s0  is  the  dielectric  permittivity  of  the  vacuum,  and  e,-  is  the  static  dielectric  constant  of  the  host 
semiconductor.  f(K)  is  a (dimensionless)  fimction,  the  exact  form  of  which  depends  on  the  details  of  the  concrete 
mechanism  leading  to  the  reduction  of  the  thermal  ionization  energy.  EAo  can  be  identified  with  the  optical  ionization 
energy  as  determined  e.  g.  from  PL  measurements. 

Several  theoretical  models  are  supporting  the  form  of  the  dependence  of  EA  on  various  parameters  especially  on 
Na  and  K as  given  by  Eq.  (1).  These  include  among  others  the  model  based  on  the  concept  of  electrostatic  interaction 
between  the  free  carriers  and  the  oppositely  charged  impurity  centres  originally  proposed  by  Pearson  and  Bardeen'®  and 
modified  by  Debye  and  Con  well,' ' a similar  model  incorporating  also  the  effects  of  the  relaxation  of  the  distribution  of 
carriers  on  the  impurity  centres  due  to  Monecke  et  al.,'^  and  another  model  based  on  the  concept  of  random  potential 
fluctuations  and  quantum  overlap  effects  changing  the  energy  spectrum  of  impurities,  proposed  by  Lien  and  Shklovskii'® 
(see  also  references'7*'").  Ajj  [eacj  £q  (j)  but  with  different  functional  forms  of  f(K).  Notably,  according  to  Pearson 
and  Bardeen'®  f(K)  = 1,  according  to  Debye  and  Con  well' ' f(K)  = 1.646(47^3)'^  = 2.6533K'^,  according  to  Monecke 
et  al.  f(K)  = F(2/3)(47t/3)E3  = 2.1828K'y,3,  and  according  to  Lien  and  Shklovskii ' ^ f(K)  is  a more  complicated  function 
with  the  value  of  1 at  K = 0,  the  value  increasing  slowly  up  to  about  1 .4  for  K = 0.4  - 0.6,  and  dropping  more  sharply  for  K 
> 0.8,  and  even  changing  sign  above  K = 0.93. 

3.  ANALYSIS  OF  THE  HALL  DATA  ON  p-GaN:Mg  AND  DISCUSSION 

Using  the  Hall  data  for  GaN:Mg  collected  from  references  '9-28  j calculated  f(K)eXp  i.e.  the  "experimental"  values  of  the 
function  in  Eq.  (1)  for  various  samples.  In  order  to  have  a meaningful  comparison  the  hole  concentration  versus  reciprocal 
temperature  curves  were  reanalyzed  using  the  standard  partially  compensated  single  acceptor  model: 

P(P  + ND)  1 ea 

= — Nv  exp( ) (2) 

NA-NA-p  p kT 

with  Nv  = 4.829xl0'5(m*/mo)3/,2T3/2  The  adjustable  parameters  were  the  acceptor  concentration  NA,  the  donor 
concentration  Njy,  and  the  thermal  ionization  energy  of  the  acceptors  EA.  The  values  used  for  the  constants  were  as 
follows:  for  the  impurity  level  spin  degeneracy  B,  a value  of  4 was  chosen,  and  the  density-of-states  effective  mass  of  holes 
in  the  valence  band  was  assumed  to  be  m /m0  = 0.8. 


Compensation  degree  K = N0/NA 

Fig.  1.  Theoretical  and  "experimental"  curves  of  the  function  f(K)  versus  K.  Data  collected  and  reanalyzed  from  the 
literature. ' ®‘~7  Theoretical  curves:  A - according  to  Debye  and  Conwell,"  B - according  to  Monecke  et  a).,'7  C - 
according  to  Lien  and  Shklovskii. 
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Fig.  2.  Plot  of  the  values  of  a versus  the  reciprocal  dielectric  constant  for  various  semiconductors.  The  straight  line 
corresponds  to  as,.  = 38.2xl0'5  meVcm.  Data  and  plot  from  Podor,13’14  the  GaN  datum  point  is  the  present  result. 


The  results  obtained  for  f(K)  versus  K are  presented  in  Fig.  1 together  with  the  theoretical  curves  for  the  models  of 
Debye  and  Conwell,11  Monecke  at  al.,12  and  Lien  and  Shklovskii.16  In  evaluating  f(K)  from  Eq.  (1)  the  acceptor 
ionization  energy  at  infinite  dilution  was  taken  as  EAo  = 224  meV,  and  the  relative  dielectric  constant  as  Sj.  — 9.6. 
Notwithstanding  the  relatively  big  scatter,  the  data  follow  a common  trend,  i.  e.  the  values  of  f(K)exp  deduced  from  the 
Hall  curves  reported  in  the  literature  increase  monotonously  with  increasing  K,  and  the  actual  values  of  the  function 
f(K)exn  seem  to  be  somewhat  better  described  by  the  model  due  to  Debye  and  Conwell  1 1 than  by  the  model  of  Monecke  et 
al.,  ^without  any  fitting  parameter.  It  is  to  be  noted  that  in  the  limit  of  K = 0 the  deduced  values  of  f(K)exp  tend  to  zero. 
The  theoretical  curve  based  on  the  model  of  Lien  and  Shklovskii^  is  clearly  at  variance  with  the  data  analysed  here.  In 
the  range  of  intermediate  compensation  (K  = 0.3  to  0.6)  the  values  of  f(K)exp  lie  significantly  higher  than  the  curve 
corresponding  to  this  model,  and  at  high  compensation  degrees  this  difference  reaches  a factor  of  2 to  3.  Similar  behaviour 
of  the  function  f(K)  in  various  semiconductor-impurity  systems  has  already  been  noted  by  Zabrodskii  and  Timofeev1  ^ and 
by  the  present  author15  too. 

If  we  suppose  the  a priori  applicability  of  the  electrostatic  interaction  model  of  Debye  and  Conwell  to  the 
experimental  data,  as  the  validity  of  this  assumption  is  being  corroborated  by  the  results  of  analysis  of  a great  amount  of 
data,1®-14  then  using  the  relationship  Nq  = KNA  Eq.  (1)  can  be  rewritten  as 


e2 

EA  = EAo  - 2.6533 ND1/3  = EAo  - aND1/3  (3) 

47terso 


with  a theoretical  value  of  a = 38.2x10  5/ej.  (meVcm).  Using  the  model  of  Monecke  et  al.12  the  constant  would  be  a - 
3 1 .4x10"5/£,.  (meVcm). 
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Fig.  3.  Plot  of  the  thermal  activation  energy  of  Mg  acceptors  in  GaN  versus  the  cubic  root  of  the  compensating  impurity 
(donor)  concentration.  Up  triangles,  solid  circles,  down  triangles,  squares,  circles,  solid  squares  - data  from  the 
literature.  19,21,22,26,27,28  The  straight  line  represents  the  least-squares  fit  to  the  data. 

Fig.4.  Acceptor  activation  energies  deduced  using  a model  incorporating  the  dependence  of  the  activation  energy  on  the 
concentration  of  ionized  acceptors:  Naj  = Np  + p.  Triangles  - data  from  Nakayama  et  al.,21  squares  - data  from  Gotz  et 
al.,23  reanalyzed  with  EAo  = 224  meV. 

This  model  suggests  that  the  value  of  a for  various  impurities  in  the  same  semiconductor  host  should  be  equal,  and 
the  values  of  a for  different  semiconductors  should  scale  with  the  reciprocal  of  the  dielectric  constants.  An  analysis  of  the 
available  large  amount  of  data  for  group  IV,  III-V,  II- VI  and  II-IV-V2  semiconductors  confirms  this  statement.13’14  Fig. 
2,  after13’14  summarizes  these  results.  Empirically  it  was  found  that  a = (40±5)x10"5/Sj.  meVcm,  which  corresponds  to 
the  theoretical  value  for  the  model  of  Debye  and  Conwell.1 1 Fig.  3.  presents  the  data  for  GaN:Mg  plotted  according  to  Eq. 
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(2).  Least-squares  fit  yields  EA  = (217  ± 8)  - (3.14  ± 0.40)xl0‘5ND1/3  (energy  in  meV,  concentration  in  cm'3).  The 
energy  intercept  agrees  within  error  with  the  optical  ionization  energy,  the  fitted  value  of  the  slope  a however,  as  can  be 
seen  in  Fig.  2,  is  lower  than  the  theoretical  value  for  the  model  of  Debye  and  Con  well,  ^ ^ and  is  close  to  that  expected  for 
the  model  of  Monecke  et  al.  *3  However,  if  the  fit  is  performed  with  a fixed  intercept  EAo  = 224  meV,  a greater  value  of  a 
is  obtained,  which  is  close  to  the  value  expected  on  the  basis  of  the  model  of  Debye  and  Conwell.  ^ * 

The  usual  interpretation  attached  to  Eq.  (3)  is  that  the  thermal  ionization  energy  decreases  linearly  with  the  cubic 
root  of  the  ionized  acceptor  (majority  impurities  in  p-type  material)  concentration  because,  at  least  at  low  temperatures 
where  p «NA,  Npj,  NAj  ~ Nj).  However,  in  general  NA;  = Np  + p,  therefore  in  the  strict  sense  of  the  electrostatic 
interaction  models  (like  of  Debye  and  Conwell  “and  of  Monecke  et  al.^3)  in  Eq.  (3)  NAj  = Nq  + p should  figure  instead 
of  With  this  refinement  instead  of  Eq.  (3)  we  have 

Ea  = EAo  - aNAi,/3  = EAo  - a(ND  + p)1/3  (4) 

implying  that  the  "effective"  thermal  ionization  energy  itself  depends  on  the  temperature  through  the  temperature 
dependence  of  the  carrier  concentration  p.  Typical  results  of  Hall  curve  fits,  taking  into  account  the  implicit  temperature 
dependence  of  the  thermal  activation  energy  as  given  by  Eq.  (4),  are  presented  in  Fig.  4.  In  these  fits  EAo  =224  meV,  a 
constant,  was  assumed.  As  might  be  expected  from  Eq.  (4),  the  fitted  activation  energies  exhibit  a substantial  decrease  at 
higher  temperatures,  and  concomitantly  and  somewhat  surprisingly,  the  fitted  values  of  the  constant  a are  20-40  per  cent 
higher  than  those  obtained  from  the  conventional  analysis,  and  also  from  the  electrostatic  interaction  model.  ^ 


4.  CONCLUSIONS 

The  observed  differences  in  the  thermal  and  optical  ionization  energies  of  Mg  acceptors  in  GaN  were  analyzed  on  the  basis 
of  the  available  models  for  the  concentration  and  compensation  degree  dependence  of  thermal  ionization  energies  of 
impurity  centres.  It  was  established  that  these  differences  can  be  satisfactorily  described  using  a simple  electrostatic 
interaction  model  originally  due  to  Pearson  and  Bardeen  and  Debye  and  Conwell. 


ACKNOWLEDGEMENTS 

This  work  was  supported  from  grant  No.  030395  from  the  Hungarian  National  Research  Fund  (OTKA). 


REFERENCES 

1.  S.  Strite  and  H.  Morkoc,  "GaN,  AIN,  and  InN:  a review,"  J.  Vac.  Sci.  Technol.  B 10,  pp.  1237-1266,  1992. 

2.  S.  Strite,  M.  E.  Lin  and  H.  Morkoc,  "Progess  and  prospects  for  GaN  and  the  III-V  nitride  semiconductors,"  Thin  Solid 

Films  231,  pp.  197-210,  1993. 

3.  J.  W.  Orton  and  C.  T.  Foxon,  "Group  Ill-nitride  semiconductors  for  short  wavelength  light  emitting  diodes,"  Rep. 

Prog.  Phys.  61,  pp.  1-75,  1998. 

4.  S.  J.  Pearton,  J.  C.  Zolper,  R.  J.  Shul  and  F.  Ren,  "GaN:  Processing,  defetcs,  and  devices,"  J.  Appl.  Phys.  86,  pp.  1-78, 

1999. 

5.  H.  Amano,  H.  Kito,  K.  Hiramitsu  and  I.  Akasaki,  "P-type  conduction  in  Mg-doped  GaN  treated  with  low-energy 

elecrtron  beam  irradiation,"  Jpn.  J.  Appl.  Phys.  28,  L21 12-L21 14, 1989. 

6.  S.  Nakamura,  T.  Mukai,  M.  Senoh  and  N.  Iwasa,  "Thermal  annealing  effects  on  p-type  Mg-doped  GaN  films,"  Jpn.  J. 

Appl.  Phys.  31,  pp.  L139-L142,  1992. 

7.  M.  A.  L.  Johnson,  Zh.  Yu,  C.  Boney,  W.  C.  Hughes,  J.  W.  Cook,  J.  F.  Schetzina,  H.  Zhao,  B.  J.  III-V  Nitrides, 

Materials  Research  Society  Symposium  Proceedings,  edited  by  F.  A.  Ponce,  T.  D.  Moustakas,  I.  Akasaki  and  B.  A. 
Monemar,  Vol.  449,  pp.  215-220,  MRS,  Pittsburgh,  1997. 

8. 1.  Akasaki,  K Amano,  M.  Kito  and  K.  Hiramitsu,  "Photoluminescence  of  Mg-doped  p-type  GaN  and 

electroluminescence  of  GaN  p-n  junction  LED,"  J.  Lumin.  48/49,  pp.  666-670,  1991 . 

9. 1.  Akasaki  and  H Amano,  "MOVPE  growth  of  high  quality  AlxGaj_xN/GayInj_yN  heterostructures  for  short 

wavelength  light  emitter,"  in  Materials  Research  Society  Symposium  Proceedings,  Vol.  339,  pp.  443-451,  MRS, 
Pittsburgh,  1994. 


302 


Proc.  SPIE  Vol.  4412 


10.  G.  L.  Pearson  and  J.  Bardeen,  "Electrical  properties  of  the  silicon  and  silicon  alloys  containing  boron  and  phosphorus," 

Phys.  Rev.  75,  pp.  865-883,  1949. 

1 1.  P.  P.  Debye  and  E.  M.  Conwell,  "Electrical  properties  of  n-type  germanium,"  Phys.  Rev.  93,  pp.  693-706,  1954. 

12.  J.  Monecke,  W.  Siegel,  E.  Ziegler  and  G.  Ktihnel,  "On  the  concentration  dependence  of  the  thermal  impurity-to-band 

activation  energies  in  semiconductors,"  phys.  stat.  sol.  (b)  103,  pp.  269-279,  1981. 

13.  B.  Podor,  "On  the  concentration  dependence  of  the  thermal  ionization  energy  of  impurities  in  InP,”  Semicond.  Sci. 

Technol.  2,  pp.  177-178,  1987. 

14.  B.  Podor,  "On  the  concentration  dependence  of  the  thermal  ionization  energy  of  impurities  in  semiconductors," 

Proceedings  of  the  2nd  International  Institute  on  New  Developments  in  Semiconductor  Physics,  August  30  - 
September  4,  1987,  Szeged,  Lecture  Notes  in  Physics,  edited  by  G.  Ferenczi  and  F.  Beleznay,  Vol.  301,  pp.  55-60, 
Springer  Verlag,  Berlin,  1988. 

15.  B.  Podor,  "On  the  concentration  and  composition  dependence  of  the  acceptor  energy  in  AlxGaj.xAs:Ge  (x<0.40)," 

phys.  stat.  sol.  (a)  119,  pp.  K135-K139,  1979. 

16.  N.  V.  Lien  and  B.  I.  Shklovskii,  "Level  of  percolation  in  weakly  doped  semiconductors,"  Fiz.  Tekh.  Poluprov.  13,  pp. 

1763-1770, 1979. 

17.  B.  I.  Shklovskii  and  A.  L.  Efros,  Electronic  Properties  of  Doped  Semiconductors,  Springer-Verlag,  Berlin,  1984. 

18.  A.  G.  Zabrodskii  and  M.  P.  Timofeev,  "On  the  effect  of  random  field  on  thermal  energy  of  ionization  in  lightly  doped 

semiconductors,"  Fiz.  Tekh.  Poluprov.  21,  pp.  2217-2219, 1987. 

19.  T.  Tanaka,  A.  Watanabe,  H.  Amano,  Y.  Kobayashi,  I.  Akasaki,  S.  Yamazaki,  and  M.  Koike,  "P-type  conduction  in 

Mg-doped  GaN  and  A1q  QgGaQ  grown  by  metalorganic  vapor  phase  epitaxy,"  Appl.  Phys.  Lett.  65,  pp.  593- 
594,1994. 

20.  W.  Gotz,  N.  M.  Johnson,  J.  Walker,  D.  P.  Bour,  FL  Amano  and  I.  Akasaki,  "Hydrogen  passivation  of  Mg  acceptors  in 

GaN  grown  by  metalorganic  chemical  vapor  deposition,"  X/?/?/.  Phys.  Lett.  67,  pp.  2666-2668,  1995. 

21.  H.  Nakayama,  P.  Hacke,  M.  R.  H.  Khan,  T.  Detchprohm,  K.  Hiramitsu  and  N.  Sawaki,  "Electrical  transport  properties 

of  p-GaN,"  Jpn.  J.  Appl.  Phys.  35,  pp.  L282-L284,  1996. 

22.  W.  Gotz,  N.  M.  Johnson,  J.  Walker,  D.  P.  Bour  and  R.  A.  Street,  "Activation  of  acceptors  in  Mg-doped  GaN  grown  by 

metalorganic  chemical  vapor  deposition,"  Appl.  Phys.  Lett.  68,  667-669, 1996. 

23.  W.  G8tz,  N.  M.  Johnson,  D.  P.  Bour,  C.  Chen,  FL  Lin,  C.  Kuo  and  W.  Umer,  "Shallow  and  deep  level  defects  in  GaN," 

in  Gallium  Nitride  and  Related  Materials,  Materials  Research  Society  Symposium  Proceedings,  edited  by  F.  A. 
Ponce,  R.  D.  Dupuis,  S.  Nakamura  and  J.  A.  Edmond,  Vol.  395,  pp.  443-454,  MRS,  Pittsburgh,  1996. 

24.  Y.  Ohuchi,  K.  Tadatomo,  FL  Nakayama,  N.  Kaneda,  T.  Detchprohm,  K.  Hiramitsu  and  N.  Sawaki,  "New  dopant 

precursors  for  n-type  and  p-type  GaN,"  J.  Cryst.  Growth  170,  pp.  325-328,  1997. 

25.  W.  Kim,  A.  E.  Botchkarev,  A.  Salvador,  G.  Popovici,  H.  Tang  and  H.  Morkoc,  "On  the  incorporation  of  Mg  and  the 

role  of  oxygen,  silicon,  and  hydrogen  in  GaN  prepared  by  reactive  molecular  beam  epitaxy,"  J.  Appl.  Phys.  82,  pp. 
219-226, 1997. 

26.  B.  Suchanek,  M.  Palczewska,  M.  Pakula,  J.  Baranowski  and  M.  Kaminska,  "Electrical  and  ESR  studies  of  GaN  layers 

grown  by  metal  organic  chemical  vapour  deposition,"  Acta  Phys.  Polonica  A 92,  pp.  1001-1004,  1997. 

27.  M.  S.  Brandt,  P.  Herbst,  FL  Angerer,  O.  Ambacher,  and  M.  Stutzmann,  "Thermopower  investigation  on  n-  and  p-type 

GaN,"  Phys.  Rev.  B 58,  7786-7791,  1998. 

28.  P.  Kozody,  H.  Xing,  S.  P.  DenBaars,  U.  K.  Mishra,  A.  Sa?der,  R.  Perrin,  S.  Elhamri  and  W.  C.  Mitchell,  "Heavy 

doping  effects  in  Mg-doped  GaN,"  J.  Appl.  Phys.  87,  pp.  1832-1835, 2000. 

29.  W.  Siegel,  G.  Kfihnel,  H.  Koi  and  W.  Gerlach,  "Electrical  properties  of  n-type  and  p-type  InP  grown  by  the  synthesis, 

solute  diffusion  technique,"  phys.  stat.  sol  (a)  95,  pp.  309-316,  1986. 


Proc.  SPIE  Vol.  4412 


303 


